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1. Introduction 

Top-quark production is one of the most important processes at hadron colliders. Within 
the Standard Model, top quarks can be produced in pairs (via strong interaction) or in- 
dividually (via electroweak processes). Top-antitop pair production is known to have the 
largest cross section, and has been studied extensively, both experimentally and theoreti- 
cally. In fact, the top quark was discovered in 1995 at the Tevatron in ti production events. 
The observation of single-top electroweak production is instead more difficult and it was 
announced by the Tevatron experiments more recently [1,2]. This was due not only to the 
fact that single-top cross section is smaller than the ti one, but also to the presence of large 
backgrounds, namely W + jet and ti, that required highly non-trivial analysis strategies. 

Although experimentally challenging, single-top electroweak production is particularly 
important since it provides a relatively clean place to study the electroweak properties 
of the top quark. For instance, it allows a direct measurement of the Vtb CKM matrix 
element [3], which is important in testing the unitarity of the CKM matrix. Moreover, 
electroweak-produced top-quarks are highly polarized, hence angular correlations of top- 
quark decay products are particularly sizeable [4,5], providing a good probe of the spin 
structure of the W-t-b vertex. Finally, since the top-quark mass is close to the scale of 
electroweak symmetry breaking, new physics effects could be discovered in the the top 
quark sector, and reactions where only a single-top is produced are particularly sensitive 
to some BSM models (see, for example, refs. [6,7]). Single-top production (in the Wt- 
channel) is also an important background for some Higgs-boson search channels, such as 
H -»• W+W [8]. 

The hadroproduction of single-top quarks in proton-proton collisions is traditionally 
classified according to the partonic processes present at LO: the s-channel processes (qq — >• 
tb) involve the exchange of a time-like W boson, the i-channel processes (bq — > tq') involve 
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the exchange of a space- like W boson, while associated Wt production (bg — > tW~) involves 
the production of a top quark in association with a W boson. The t-channel process is 
the main source of single-top quarks, both at the Tevatron and the LHC. At the Tevatron 
the Wt contribution is negligible and the s-channel cross section is roughly half of the t- 
channel's one. At the LHC, instead, the VFi-production cross section is a factor 3 less than 
the i-channel's one, while the s-channel is negligible. Therefore, an accurate description of 
all the three production channels is important: in particular, at the LHC, the H^t-channel 
will play a significant role. 

Given the above reasons, it is then desirable to reach a high precision in the theoretical 
predictions, both for total rates and for more exclusive distributions. Higher order correc- 
tions to single-top hadroproduction have been calculated in refs. [9-20]. Nowadays, one of 
the ways to go beyond this level of accuracy is to merge the fixed order accuracy of a NLO 
calculation with the (Next-to)-Leading-Logarithm accuracy (and the flexibility) of a (par- 
ton shower) Monte Carlo event generator. At present there are two methods to interface 
NLO calculations with parton showers in a consistent way: MCONL0 [21] and POWHEG [22,23]. x 
They have been applied to several processes at lepton [27-29] and hadron [30-44] collid- 
ers, in conjunction with the HERWIG [45], PYTHIA [46] and HERWIG++ [47] parton shower 
algorithms. 2 3 Despite the theoretical formulation of the two methods being quite differ- 
ent, reasonable agreement has been found as well as the reason for the (few) differences 
encountered (see for example [50] for a recent discussion). 

In ref. [41] the details of the implementation of single-top s- and t-channel production 
processes in the POWHEG framework have already been described. The aim of this paper is 
to describe the implementation of the V7t-channel production process. Therefore, results 
presented here can be seen as a completion of the work described in the previous reference. 

This paper is organized as follows. In sec. 2 we describe some of the technicalities of the 
implementation. We also summarize how the problem of the interference between the NLO 
corrections to this process and the ti calculation has been dealt with. Here we anticipate 
that we have used the same strategy adopted by the MC@NL0 authors, i.e. two definitions 
have been implemented, such that their difference can be considered a measure of one of 
the theoretical uncertainties that affects this process. In sec. 3 results are presented: in 
particular we compare POWHEG and MC@NL0 results for the two aforementioned definitions of 
the NLO corrections. We also show results obtained with typical cuts and results obtained 
using the PYTHIA shower. Finally, in sec. 4 we give our conclusions. 

2. The POWHEG implementation 

To match single-top VFt-channel production NLO corrections with a parton shower using 
the POWHEG method, the P0WHEG-B0X package has been used. P0WHEG-B0X is a program 

1 In literature, other proposals exist [24-26]. The MCONLO and POWHEG methods are the only ones where, 
currently, several full processes at lepton and hadron colliders have been implemented. 

2 In [48] a compared study of the effects of different parton shower algorithms has also been performed, 
in case of Higgs-boson production via gluon fusion. 

3 An implementation of the POWHEG method has also appeared recently within the SHERPA event genera- 
tor [49]. 
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that automates all the steps described in ref. [23], turning a NLO calculation into a PDWHEG 
simulation. The details of how the program works have been largely described in ref. [51]. 
In practice, starting from some of the typical building blocks of the NLO calculation for 
the process at hand, the program produces a set of partonic events ready to be showered 
by a shower Monte Carlo program. 

In this section we present some of the inputs that we calculated and then coded in the 
format needed by the POWHEG-BOX package to work. In particular, we describe how the LO 
kinematics and the relevant differential cross sections up to next-to-leading order in the 
strong coupling a s have been obtained. All quark masses have been set to zero (except, 
of course, the top-quark mass) and the full Cabibbo-Kobayashi-Maskawa (CKM) matrix 
has been taken into account in the calculation. However, for sake of illustration, through 
this section we set the CKM matrix equal to the identity. Therefore, the d-type quark 
connected to the top and the W boson will be denoted as b. Furthermore, in this paper 
we always refer to top-quark production: anti-top production is obtained simply by charge 
conjugation. 

With this convention, the LO partonic process for Wt production is 



We denote with Bf, g the (summed and averaged) squared amplitude, divided by the flux 
factor. The process initiated by the partons gb is also present. For brevity, throughout 
the paper processes that can be obtained from the written ones by simply exchanging the 
order of the incoming partons will be omitted from the formulae. 

Before giving the relevant formulae, here we stress again that the NLO corrections 
to the Wt production channel are not well defined, due to interference effects with the ti 
process. This problem is well known, and several approaches have been introduced [15, 
52,53]. To deal with this theoretical issue, we used the same strategy first described 
in the corresponding MCSNLO publication [36] and later extensively studied in ref. [54]: 
two definitions for the NLO corrections have been considered, and the relative difference 
between the two can be interpreted as a measure of the theoretical uncertainty in the 
definition of genuine NLO corrections to the VFi-channel process. Since this issue concerns 
the radiative part of the QCD NLO corrections, a more accurate discussion will be given 
in sec. 2.3. 

2.1 Born kinematics 

Following the notation of ref. [23], we denote with k 9 and k e the incoming parton momenta, 
aligned along the plus and minus direction of the z axis, and by k\ and k,2 the outgoing 
IF-boson and top-quark momenta, respectively. The top-quark and VF-boson masses are 
denoted by mt and mw- During the POWHEG-algorithm step where the inclusive NLO cross 
section is evaluated (the calculation of the B function), their virtualities are kept fixed and 
equal to their masses. If K 9 and K e are the momenta of the incoming hadrons, then we 
have 



b + g ->• W~+t. 



(2.1) 




(2.2) 
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where the momentum fractions, and momentum conservation reads 

k 9 + k e = h + k 2 . (2.3) 

We now introduce the variables 

'=<^ H ^ttS'-SlS - (2 - 4) 

and 8, the angle between the outgoing top quark and the fc e momentum, as seen in the 
partonic center-of-mass (CM) frame. We denote with 4> the azimuthal angle of the outgoing 
top quark in the same reference frame. Since the differential cross sections do not depend 
on the overall azimuthal orientation of the outgoing partons, we set this angle to zero. At 
the end of the generation of an event, we perform a uniform, random azimuthal rotation 
of the whole event, in order to cover the whole final-state phase space. The set of variables 
$2 = {s, Y, 8, 0} fully parameterizes the Born kinematics. From them, we can reconstruct 
the momentum fractions 

* Y - I & -Y 



x ® = \l~s e } x ° = y~s e ' ^ 2 ' 5 ^ 

where S = (K e + K e ) 2 is the squared CM energy of the hadronic collider. The outgoing 
momenta are first reconstructed in their longitudinal rest frame, where Y = 0. In this 
frame, their energies are 

-k° 1 W-o= ~ S + r f- m2t and k%-_ = ~ S + m y r m ™ . (2.6) 
2\/s 2ys 

The two spatial momenta are opposite and their modulus \k\y=o ls found by using the 
on-shell condition = (fc°|y =0 ) — (|fc|y =0 ) 2 . We fix the top-quark momentum to form 
an angle 8 with the © direction and to have zero azimuth (i.e. it lies in the xz plane and 
has positive x component). Both k\ and k-2 are then boosted back in the laboratory frame, 
with boost rapidity Y. The Born phase space, in terms of these variables, can be written 
as 

d* 2 = dx, dx 6 (2^(k. + k 6 -h- h) (^j||o 

= i— ds dY dcos8 ^, (2.7) 
S 16vr 2vr y J 



where 



^=VW, with 2(m^ + m?) _ (m^-r^)^ 

s s z 



The content of eq. (2.7) and the procedure described above to define the parton momenta 
from the variables set $2 have been coded in the subroutine born_phsp of the P0WHEG-B0X 
package [51]. 
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Figure 1: Feynman diagrams for the LO process bg — > W t. The top quark is denoted with a 
double line. 

2.2 Born and virtual contributions 

The squared matrix element for the lowest order subprocess bg —> W~t has been obtained 
with MadGraph [55]. It has been checked with the expression reported in eq. (3.5) of 
ref. [36] and with an independent calculation performed with FeynCalc [56], starting from 
the two Feynman diagrams in fig. 1. The top-quark width Tt has been set here to zero. 

To regularize soft and collinear divergences, the POWHEG-BOX package uses an auto- 
mated implementation of the subtraction algorithm proposed by Frixione, Kunszt and 
Signer [57,58]. The counterterms needed for the numerical subtraction of the real squared 
amplitudes soft divergences are obtained from the color-linked Born squared amplitudes 
Bij defined in eq. (2.97) of ref. [23], using the eikonal approximation. For the process at 
hand, the By matrix is proportional to the full Born squared amplitude. In fact, since at 
the lowest order we have only three color-connected partons, the equality 

T e + T e + T 2 = (2.9) 

holds, where Tj is the color operator associated with the i parton in the Born process (2.1). 
Using the property 

TiTj J T i + Tj )^T ? -T^ (2io) 

and the fact that T| = C F if i is a quark and T| = C A if i is a gluon, we obtain Bij = 
(2C F — C A )Bb g /2 if neither i or j are gluons, and By = C A Bb g /2 if i or j is a gluon. 

Collinear counterterms are obtained using the collinear factorization. However, since in 
the VFt-channel LO process (2.1) an external gluon is present, the collinear limits associated 
with this leg do not factorize in terms of the Altarelli-Parisi unpolarized splitting kernels 
times the Born contribution Bb g - In fact azimuthal correlations in the branching process 
are present, and to build a local counterterm the POWHEG-BOX package makes use of the spin 
correlated Born cross sections B^ IU , defined in eq. (2.8) of ref. [23]. The FeynCalc program 
has been used to calculate this matrix and to translate the result in a Fortran routine. 

One loop virtual contributions have been calculated and algebraically reduced to scalar 
integrals using the Passarino-Veltman algorithm. We used the same renormalization pro- 
cedure described in ref. [36]. We checked that the IR pole structure we are left with 
corresponds to the singularities of the real contributions. The finite part that enters in 
the soft-virtual contribution Vb g (eq. (4.6) of ref. [51]) is computed numerically, using the 
package QCDloop [59] to evaluate the finite part of the scalar integrals. 
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2.3 Real contributions 



The real emission corrections can be classified as follows: 



b + q(q) ->■ W 
g + g -> W 
g + <? -> IF 



6 + 5 ->■ IF 



+ t + q(q) 
+ t + b, 



(2.11) 
(2.12) 
(2.13) 
(2.14) 



We denote the corresponding contributions to the cross section as TZb g , lZb q , 1Z gg and lZ qq - 
respectively, where we used again the standard POWHEG notation first introduced in ref. [23]. 
The 3-body phase space is denoted as <&3 and the corresponding momenta as k s , k e , k\, 
ki and £3, where ^3 is the momentum of the outgoing light parton (the FKS parton) while 
the other momenta correspond to those of the two incoming partons, the VF-boson and the 
top-quark. In fig. 2 some representative diagrams are shown. 



Figure 2: Representative Feynman diagrams for the processes (2.11), (2.12) (up), (2.13) (center) 
and (2.14) (bottom). 

The processes (2.13) and (2.14) have a final state that corresponds to ti production 
followed by the decay i —> W~b. A consequence of this fact is the well known problem of 
interference between Wt and ti production and it is the reason why QCD NLO corrections 
to the VFi-channel are not well defined. In the following we will state more precisely the 
nature of the problem and explain how we dealt with it. 

As one can see from the last four Feynman diagrams in fig. 2, diagrams associated with 
the subprocesses (2.13) and (2.14) can be divided into two sets. Following the nomenclature 
of ref. [36], we call doubly-resonant the diagrams where a top-quark propagator goes on- 
shell when the invariant mass of the system made by the IF-boson and the outgoing 6-quark 
( m Wb) approaches nif. The other diagrams we call singly-resonant. Since these two sets 
of diagrams have to be summed at the amplitude level, interference effects are present. 
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This interference between NLO real corrections to Wt and lowest order ti production 
(followed by a decay) would not be a problem if the contamination on tt-like diagrams 
was numerically negligible with respect to the size of singly-resonant diagrams. This is 
certainly not the case when one approaches the region (kw + ki) 2 = (k\ + k^) 2 ~ m 2 . 
In fact, in this region the real corrections to doubly-resonant diagrams become huge (an 
internal propagator goes on-shell), and the perturbative expansion (in power of a s ) for the 
NLO corrections to single-top W^i-channel loses its meaning. 
To deal with this problem, several approaches are possible: 

• The more drastic approach is to consider that top quarks are not detectable parti- 
cles. This approach would remove any interference issue, since processes would be 
unambiguously classified accordingly only to experimentally measurable (QCD) final 
states: one would have the processes W + W~bb and W~W + b. At present, the price 
to pay would be to neglect NLO corrections, since these are not available for processes 
where top-quarks are not on-shell. 4 

• An alternative approach is to give a prescription for removing the contamination from 
it-like contributions, while keeping the top-quark as a final state particle. This task 
can be accomplished using cuts to avoid the doubly-resonant region or including an 
ad hoc subtraction term to remove (or suppress) the ti contribution. 

The cut strategy was first adopted in ref. [52], where an explicit cut on m w i was 
used. In ref. [15] a 6-jet veto was instead used, together with a careful choice of the 
factorization scale. Moreover, in ref. [15] for the first time differential results at the 
NLO were presented and QCD corrections to the top-quark decay have been also 
included, in the narrow-width approximation. 

In ref. [53] a subtraction term was first introduced in the context of an inclusive 
LO(+LL) calculation. A similar method has been adopted in the fully inclusive NLO 
calculation presented in ref. [13]. 

Finally, in refs. [36,54] the interference problem was reexamined at length, and it was 
shown that a separate treatment of ti and Wt production is feasible, also when the 
NLO calculation is matched with a parton shower. 

In this work we have used the same approach described by the MC@NL0 authors. This 
relies on the observation that a meaningful definition of the Wt-ch.au.nel process (as a signal 
or a background) is possible only through cuts on final state objects. If interference effects 
with ti are negligible after these cuts are applied, then it is possible to consider VFt-channel 
a well defined process. Since cuts act differently in different phase-space regions, one needs 
to quantify the interference between Wt and ti locally in the phase-space. To this end, two 
definitions for the NLO corrections were given. It was shown that by comparing the two 

4 We also recall that a calculation of NLO corrections to the production and the decay of top quarks 
performed in the narrow width approximation (as the one appeared in ref. [15]) would not avoid the 
interference problem. To our knowledge, the only calculations that fully include top offshellness effects are 
the ones described in ref. [60, 61] for single-top s- and t-channel production respectively. 
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results it is possible to address the previous question, and it was found that, for several 
sets of cuts, the theoretical uncertainty due to interference effects is typically negligible 
with respect to other theoretical errors. In the following we will discuss briefly the two 
subtraction strategies and how they have been implemented within the POWHEG framework. 
In sec. 3 the corresponding results will be shown. 

The two definitions proposed in ref. [36] are known as Diagram Removal (DR) and 
Diagram Subtraction (DS). Their difference can be better understood by writing a generic 
amplitude for the processes (2.13) and (2.14) as 

M =M wt + M a , (2.15) 

where Ai wt and Ai u denote respectively the sum of all the singly- and doubly-resonant 
Feynman diagrams for the partonic subprocess at hand. In DR one defines the real contri- 
bution 1Z by eliminating the tt contribution M tt from M before squaring the amplitude. 
Instead in DS one keeps the full squared amplitude but subtracts from it a local countert- 
erm C SUB in order to suppress the tt contribution at the cross section level. In this respect, 
DS can be seen as a refinement of the method proposed in ref. [53]. Schematically, we have: 

\M wt \ 2 

K OR = 1 1 , (2.16) 

\M wt + M tl \ 2 - C SUB 
K DS = l — , (2.17) 

where s is the squared CM energy. Some comments are due here: 

• The DR method is not gauge invariant. This issue was studied in depth by the 
authors of ref. [36] , and it was shown that the impact of gauge dependence in the DR 
calculation is numerically negligible. 

• In the DS approach one wants to build a gauge invariant subtraction term that exactly 
cancels the tt contribution when the doubly-resonant region is approached. Thus, the 
subtraction term C SUB has to fulfill the following requirements: 

1. gauge invariance. 

2. match exactly the |.M**| 2 contribution when the doubly-resonant region is ap- 
proached. 

3. possibly fall off (quickly) far from the doubly-resonant region. 

The third requirement is needed to keep the full NLO corrections unmodified away 
from the ti peak. Apart from the three requirements above, there is some freedom 
in the definition of C SUB . 

• By taking the difference between eq. (2.17) and (2.16), one finds 

X+ \M tl \ 2 - C SUB 
^ds _ k dr = , (2.18) 
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where X = 2S(M M" ). Therefore, the difference between results obtained with 
DR and DS can be interpreted as a measure of the size of the interference X, provided 
that the difference \M tl \ 2 - C SUB is small. 

We implemented in POWHEG the two subtraction methods. The squared amplitudes 
\A4\ 2 and [.M^ 7 *] 2 have been obtained using MadGraph. The subtraction term was chosen 
as in ref. [36]: 

™> - ((tl+t2 /:g +(m , r , )2 ix""«>i 2 . 

where <& 3 is a point in the 3-body phase space obtained by reshuffling the $3 kinematics 
in order to have (k\ + ks) 2 = m 2 , i.e. an exactly doubly-resonant configuration. Hence, the 
choice of the subtraction term is the same as the one used by the MC@NL0 authors. In fact, 
in spite of the aforementioned freedom, as it was already pointed out in ref. [36] , the choice 
of the amplitude Ai tt evaluated at the point 3> 3 is unavoidable if one wants to achieve the 
exact cancellation of the doubly-resonant contribution while retaining gauge invariance, 
since A4 ~ A4 tt when {k\ + £13) 2 — > m 2 and gauge invariance is preserved only if the 
internal t propagator is on shell. The only real freedom is in the choice of the prefactor, 
and the Breit-Wigner profile seems the more natural choice if one wants the difference 
between |.M'*| 2 and C SUB to be close to zero as much as possible away from the resonance. 

The inclusion of DR in POWHEG is straightforward. The procedure we adopted for DS 
is more subtle, since 1Z DS is not positive-definite. This affects two steps of the POWHEG 
method. In the following we describe how we proceeded in our implementation. As usual, 
we use the standard notation of ref. [23] . 

• Having a real correction that is not always positive-definite increases the chances to 
have regions of $2 where the B function becomes negative. In these cases, the radia- 
tive event generated with POWHEG starting from the underlying-Born configuration $2 
is negative-weighted. We have checked that this occurrence is rare, so that the ben- 
efits of having positive- weighted events are not spoilt. Moreover, in the POWHEG-BOX 
package, a procedure (called "folded" integration) to reduce further the occurrence 
of these events is available, as explained in sec. 4.1 of ref. [51]. By using it, we have 
verified that the occurrence of negative-weighted events can be further reduced also 
for the IFt-channel case. This also implies that the function B is positive, as in the 
other POWHEG implementations. 

Quantitatively, in the event sample we generated to produce the plots shown in 
sec. 3, the fraction of negative-weighted events was 0.05. By using the aforementioned 
folded integration, this fraction can be reduced down to 0.005, becoming therefore 
completely negligible. 

• For 3-body kinematic configurations where 1Z DS is not positive, another problem can 
occur during the generation of the hardest radiation. We recall that, to generate the 
hardest radiation, the POWHEG algorithm works by finding an upper-bound for the 
ratio TZ/B, assuming that this ratio is always positive. Negative values can therefore 
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spoil the accuracy of the method. We have checked that this happens with a certain 
frequency only close to the doubly-resonant region. Since this is the kinematic region 
where the separation of ti and Wt is already particularly critical, we decided to 
explicitly avoid to generate radiative events when close to the doubly-resonant peak. 
This has been obtained using a theta function that vanishes in this region, i.e. in DS 
the following substitution is performed when 1Z/B is evaluated: 

TZ DS -»• K DS x 6 (\m wi - m t \ - «T t ) . (2.20) 

We have tested the effects of this cutoff trying values of order 1 for k, in order to avoid 
introducing effects from this parameter in phase-space regions where it is not needed. 
The outcome of this check is that, at the end of the event generation (i.e. after the 
shower and the hadronization stage), no problems caused by this cutoff are present, 
since results depend negligibly on the value of k. The only observable where some 
dependence was observed is the differential distribution of m w i close to the doubly- 
resonant peak. When values for k in the region 3 — 5 are used, the difference with 
the MC@NL0 result is minimized. The results shown in sec. 3 have been obtained with 
k = 3. 

Although we are aware that this solution may be considered unappealing, we found 
it a reasonable choice to handle with the problem of having negative values for 1Z/B 
in the POWHEG framework. We also recall that the presence of this cutoff in our 
implementation does not affect the issue of negative weights, since k is used only in 
the generation of the hardest radiation. 

We used MadGraph to calculate the term C SUB . We also notice that to avoid the 
divergence of the internal top-quark propagator when (k\ + k^) 2 — > m 2 , in DS a non- 
vanishing value for the top width is needed. 5 

3. Results 

In this section we present our results obtained after showering with HERWIG 6.510 and 
PYTHIA 6.4.22 the partonic events generated with POWHEG. We considered top production 
at the LHC, with an hadronic center-of-mass energy VS = 10 TeV. All results have been 
obtained assuming that the top-quark decays semileptonically (t — > blv) and that the 
W^-boson involved in the hard scattering decays leptonically (W~ — > iv). Branching ratios 
have been removed, so that plots are normalized to the total cross section. 

We have used the CTEQ6M [62] set for the parton distribution functions and the 
associated value of A^ = 0.226 GeV. Furthermore, as discussed in refs. [23,32], we use 
a rescaled value A MC = 1.569 A^ in the expression for a s appearing in the Sudakov form 
factors, in order to achieve next-to-leading logarithmic accuracy. 

Although the matrix-element calculation has been performed in the massless-quark 
limit, the lower cutoff in the generation of the radiation has been fixed according to the 

5 We recall that to obtain the proper cancellation of soft and collinear divergences, M wt has to be 
calculated with T t = 0. 
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mass of the emitting quark. The lower bound on the transverse momentum for the emission 
off a massless emitter (u, d, s) has been set to the value p™ m = v5A MC . We instead 
choose p™ m equal to m c or when the gluon is emitted by a charm or a bottom quark, 
respectively. We set m c = 1.55 GeV and mi, = 4.95 GeV. 

The renormalization and factorization scales have been taken equal to the transverse 
momentum of the radiated light parton during the generation of radiation, in accordance 
with the POWHEG formalism. We have also taken into account properly the heavy-flavour 
thresholds in the running of a s and in the PDF's, by changing the number of active 
flavours when the renormalization or factorization scales cross a mass threshold. In the B 
calculation, instead, /j,r and //p have been chosen equal to the top-quark mass, whose value 
has been fixed to mj = 175 GeV. In the DS approach, the amplitudes where doubly-resonant 
graphs are present and the subtraction term C SVB have been calculated with Ft = 1.7 GeV. 

To assess the validity of the approximations and the choices we made, we compare our 
results (obtained both with DR and DS) with the MCONLO outputs. 



0.30 1 ' 1 1 1 ' 1 ' 1 1 1 1 1 1 1 ' 1 ' 1 g qq 3 — I — | — I — I — I — I — | — I — I — I — I — | — I — I — I — I — | — I — I — I — I — | — r~r 




50 100 150 200 -4 -2 2 4 

Pt [GeV] 7, T 




VF> [GeV] A0 t . w 

Figure 3: Comparisons between POWHEG (interfaced to HERWIG) and MCONLD results at the LHC 
pp collider (y/S — 10 TeV), obtained with the DR prescription. NLO results are also shown in the 
lower panel. 

In fig. 3 we show a comparison between POWHEG and MC@NL0, obtained with the DR 
prescription, without applying any cut on the final state particles. POWHEG results have 
been obtained using the HERWIG parton shower, in order to minimize differences arising 
from different shower algorithms and hadronization models. In the upper panel of fig. 3 
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we show the transverse momentum of the top quark (p^) and the pseudorapidity of the 
VF-boson (r] W ) produced in the hard process (i.e. not the VL-boson present in the decay 
chain of the top quark). As expected, we found very good agreement between the two 
results, since the shape of these two distributions is due mainly to the fixed order result. 
In the lower panel, we show instead p^^ , the transverse momentum of the system made by 
the top quark and the W boson, and A<^_jy, the difference between the azimuthal angles 
of the two particles. These two quantities are significant because their shape is affected 
by Sudakov suppression effects due to the resummation performed by parton showers. To 
stress the size of these effects, we also superimposed the fixed order (NLO) prediction to 
the last two plots. We observe good agreement between showered results, and the expected 
difference with the NLO curve, where the cancellation of soft and collinear divergences takes 
place only at the edge of the distributions (at p^ = and A(fit-w = 71 respectively). 
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Figure 4: Comparisons between POWHEG (interfaced to HERWIG) and MCONLO results at the LHC 
pp collider (^/S = 10 TeV), obtained with the DS prescription. NLO results are also shown in the 
lower panel. 



In fig. 4 we show the same set of plots, obtained with the DS prescription. The same 
considerations made above are valid also for this case. We recall that the plots shown in 
figs. 3 and 4 have been obtained using the top-quark and the W-boson momenta extracted 
from the parton shower history. Therefore, these quantities are not measurable in a real 
detector. Nevertheless, since they are useful to check the implementation, we have shown 
the corresponding results. 
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We have also included the generation of the top-quark and the VF-boson decay prod- 
ucts, according to the method originally proposed in ref. [63]. This enables the generation 
of events in which spin correlation effects in the production-decay stage are taken into ac- 
count with leading-order accuracy. In figs. 5 and 6 we show a comparison between MC@NL0 
and POWHEG for some relevant leptonic distributions. We plot the transverse momenta 
of the hardest (p^) and the second hardest charged lepton in the event. We also 

show Pt i12 \ the transverse momentum of the system made by l\ and £2, and the 
difference between the azimuthal angles of the two leptons. Here, again, we found very 
good agreement in both the DR and the DS case. We recall that a quantity like A^ 1 _^ 2 
is sensitive to spin-correlation effects, as we have verified by running the code with the 
decay-generation procedure switched off and letting the shower perform isotropic decays. 
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Figure 5: Comparisons between POWHEG (interfaced to HERWIG) and MC0NL0 results at the LHC pp 
collider (y/S = 10 TeV), obtained with the DR prescription, for leptonic quantities. 



In ref. [15] the use of a 6-veto as a method to discriminate between the Wt and 
the ti processes was proposed. This idea was also reconsidered in ref. [36], where its 
applicability was studied in the context of a NLO result merged with a parton shower. 
We have performed a similar exercise using the results obtained with POWHEG and showered 
with HERWIG. The b- veto condition is defined as follows: after having sorted in decreasing p T 
all the 6-flavoured hadrons present in the final state, an event is kept if the second-hardest 
^-flavoured hadron among those with central pseudorapidity (r/ B < 2.5) has p B < p£ veto ', 
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Figure 6: Comparisons between POWHEG (interfaced to HERWIG) and MC0NL0 results at the LHC pp 
collider (^/~S = 10 TeV), obtained with the DS prescription, for leptonic quantities. 



otherwise it is discarded. 6 In our analysis procedure we chose p^~ vcto ^ = 50 GeV. 

As expected, we observe that the difference between DR and DS results is reduced 
when the 6-veto is applied: in fact, although we have not performed a detailed study 
including also uncertainties due to scale variations, we found that the ratios between the 
total NLO cross sections using DR and DS are as follows: 



(DS) 



0.95. 



(DR) 



a 



a 



(b- veto) 
(DS) 
(6- veto) 
(DR) 



0.98. 



(3.1) 



As the MC0NL0 authors already pointed out, despite the fact that the 6-veto reduces 
the difference between DR and DS total cross sections, it is not guaranteed that all the 
differential distributions are affected in the same way. To address this question, in fig. 7 
we show the effect of the 6-veto cut for two transverse-momentum spectra. In the upper 
panel we plot the results for , obtained with the DR and the DS procedures, before 
and after imposing the 6-veto. The curves in the upper-left panel are obtained without 
cuts, those on the upper-right have been obtained keeping only the events that fulfil the 
veto condition. It can be seen that the mismatch between DR and DS in the high-p T tail 



6 If all the 6-hadrons in the central pseudorapidity region have Pt < Pt vcto \ 
6-hadron has a large pseudorapidity, the event is kept. 



or if the second-hardest 
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is less sizeable when the b-veto is applied, which is the expected result. Since in ref. [36] it 
was noticed that the transverse momentum of the system made by the two hardest leptons 
turns out to be sensitive to the treatment of the doubly-resonant region, in the lower panel 
of fig. 7 we show also the predictions for p^ li2 \ Also for this quantity the effect of the 
b- veto is to reduce the differences between DR and DS, as the plot in the lower-right panel 
shows, in accordance with what has been found in ref. [36]. 
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Figure 7: Comparisons between PDWHEG (interfaced to HERWIG) results obtained with the DR and 
DS prescriptions. Plots in the right panel have been obtained with the 6-veto described in the text. 



The conclusions of this analysis can be summarized as follows. We checked that the 
b- veto reduces the differences between DR and DS. This is the expected result, since the 
b- veto was originally proposed by the authors of ref. [15] to reduce interference effects from 
tt production. In fact, by requiring a 6-veto, one reduces the number of events with two 
hard and central 6-flavoured hadrons, which is indeed one of the typical signatures of ti 
production. The fact that DR and DS total cross sections become closer when the 6-veto is 
applied is also in accordance with the interpretation of the difference between DR and DS 
being a measure of the interference between Wt and tt production. This interpretation is 
supported by the plots in fig. 7, where it is shown that the difference between DR and DS 
is reduced also for differential distributions, when the 6-veto is applied. However, the plots 
in fig. 7 show also that observables are potentially affected by the 6-veto non uniformly 
(i.e. ratios between DR and DS results can be bin-dependent). This suggests that particular 
care should be taken when one performs a full analysis where the contribution from the Wt- 
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channel process is supposed to be relevant, since the size of interference effects may depend 
on the cuts applied. As it was already observed by the MC@NL0 authors, a comparison 
between DR and DS predictions gives an estimate of the theoretical uncertainty due to 
these effects. 

We have also performed some comparisons between DR and DS results obtained with 
the PYTHIA shower. In order to maximize consistency with the POWHEG prescriptions, we 
used the p T -ordered shower. In fig. 8 we show the pseudorapidity and the transverse 
momentum of the system made by the two hardest leptons, while in fig. 9 the transverse 
momentum of the hardest non 6-flavoured jet (p^ 1 ) and for the hardest jet that contains 
a 6-flavoured hadron (p? 1 ) are shown. The plots on the left have been obtained without 
imposing any cut. Instead, the plots on the right have been obtained using cuts similar to 
the 'Wt-signal cuts" of ref. [54]. Jets have been defined according to the k T algorithm [64], 
as implemented in the FAST JET package [65], setting R = 0.7 and imposing a lower 10 GeV 
cut on jet transverse momenta. To accept an event, we required the following properties: 

• There is exactly one 6-jet with p^ > 50 GeV and {rfl < 2.5. A 6-jet is defined as a jet 
that contains at least one 6-flavoured hadron and has p^ > 25 GeV and < 2.5. 

• There are at least two light-flavoured jets with p^ > 25 GeV and \rf\ < 2.5. The 
invariant mass of the system made by the two hardest jets among these light-flavoured 
jets has to lie within 55 and 85 GeV. 

• There is one lepton with p l T > 25 GeV and \rf \ < 2.5. This lepton has also to be 
isolated with respect to the 6-jet and the two light-flavoured jets, i.e. its distance 
from the jets in the (rj, 4>) plane has to be larger than 0.4. 

• The missing transverse energy is larger than 25 GeV. 

Although we have not performed as detailed a study as the one of ref. [54], from figs. 8 
and 9 we observe that the DR and the DS predictions are consistent (within the statistical 
accuracy) also when the above cuts are applied, as was observed in the aforementioned 
work. 
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4. Conclusions 

In this paper we have described the implementation of Wt-ch&nnel single-top production 
at next-to-leading order in QCD, in the PDWHEG framework. We have used the POWHEG-BOX 
package, which is a program that automates the algorithm first proposed in ref. [22] and 
then carefully described in ref. [23]. Having used this package, we described how we calcu- 
lated the main needed ingredients. 

Since NLO corrections to single-top production in the Wt-channel are known not to be 
well-defined (real contributions interfere with ti production) , we decided to follow the same 
strategy originally proposed by the MCONLO authors in ref. [36]: we included two definitions 
of the NLO corrections, known as DR (Diagram Removal) and DS (Diagram Subtraction), 
both of which can be used when the merging of the fixed order result with parton showers 
is performed. Moreover, the difference of results obtained with these two prescriptions 
gives an estimate of the size of interference effects. We have described how we included 
the two prescriptions in POWHEG, and how we dealt with DS and the problem of its exact 
implementation within the POWHEG method. 

To check the correctness of the whole implementation and to assess the validity of 
the choices we made, results have been compared with the MCONLO program for the LHC, 
where Wt production is relevant: we found very good agreement, both for DR and DS. We 
have also compared DR and DS results when a 6-veto is imposed. We found that a 6-veto 
reduces interference effects with ti, the difference between DR and DS results becoming 
smaller when the veto is in place. Moreover, we have also presented some results obtained 
with the PYTHIA shower and with typical defining cuts for the Wt signal. Although we have 
not performed as full an analysis as the one reported in ref. [54], good agreement between 
DR and DS has been found also when typical Wt-channel cuts are applied. 

The main purpose of this work was the completion of the work presented in ref. [41], 
where the POWHEG implementation of single-top s- and t-channel was described. Therefore, 
at present all single-top processes can be simulated in the context of a NLO+PS approach 
with POWHEG as well as with MC@NL0. We also stress that, in POWHEG, Wt-channel was the 
missing process among the ones relevant for H — > WW searches: Higgs production via 
gluon [37] and vector-boson fusion [42] are already available, as well as ti [33] and VV [66] 
production, which (together with single-top Wt) are the main backgrounds . 

The computer code for this POWHEG implementation will soon be available within the 
public branch of the POWHEG-BOX package, that can be downloaded at the site 

http : //virgilio .mib . infn. it/~nason/POWHEG 
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